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SUMMARY 


An  analytic  investigation  is  made  of  the  aerodynamic  coefficients  of 
various  "bodies  located  in  a  free-molecule— f low  field.  These  bodies  are 
the  following:  flat  plate^  cylinder,  sphere,  and  cone.  Calculations  are 
performed  using  values  of  molecular  speed  ratio  (ratio  of  stream  speed  to 
most  probable  molecular  speed)  ranging  from  0  to  20. 

The  aerodynamic  coefficients  of  a  cone  are  calculated  for  angles  of 
attack  ranging  from  0°  to  60^,  The  semi vertex  angles  of  the  cones  investi¬ 
gated  vary  from  2.5^  to  30^* 

The  calculations  are  performed  assuming  two  types  of  molecular 
reflection,  specular  and  diffuse;  in  addition,  for  the  cone,  a  third  type 
of  molecular  reflection,  wherein  impinging  molecules  are  not  re— emitted 
.from  the  body  but  are  swept  along  its  surface,  is  postulated  in  order  to 
compare  the  drag  coefficients  calculated  by  free-molecule— f low  theory  with 
values  obtained  from  continuum  theory. 


INTEODUCTION 


Consideration  of  the  problem  of  the  flight  of  high-speed  long-range 
aircraft  has  indicated  that  skin  temperature  and  drag  forces  may  be 
considerably  reduced  by  operation  at  high  altitudes.  Consequently,  it 
becomes  of  interest  to  determine  the  aerodynamic  forces  arising  from  high¬ 
speed  flight  at  high  altitudes.  Flight  at  altitudes  where  the  molecular- 
mean— free  path  is  larger  than  a  characteristic  body  dimension  is  in  the 
free-molecule— flow  regime.  Molecular-mean— free  path  as  a  function  of 
altitude  is  given  in  reference  1.  (At  an  altitude  of  75  miles,  the  mean- 
free  path  is  about  1  foot.)  Calculations  relating  skin  temperature, 
altitude,  and  velocity  have  previously  been  completed  in  reference  1; 
drag  forces,  however,  were  not  computed. 

Analytic  Investigations  of  the  drag  forces  acting  on  simple  body 
shapes  in  fVee-molecule  flow  have  been  reported  in  references  2,  3,  and 
4,  while  experimental  investigations  have  been  described  in  references 
4  and  5* 
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In  all  these  analyses^  it  vas  assumed  that  the  gas  molecules  have  a 
Maxwellian  distrihution  of  thermal  velocity  superimposed  upon  the  mass 
velocity.  This  assumption  is  likewise  made  in  the  present  paper, 

Heineman  (reference  2)  and  Ashley  (reference  3)  have  calculated  the  aero¬ 
dynamic  coefficients  for  various  body  shapes  including  the  flat  plate, 
sphere,  cylinder,  and  cone.  However,  the  expressions  for  the  drag  coef¬ 
ficients  of  the  sphere  and  cylinder  were  not  integrated  hy  Heineman, 
while  Ashley* s  solution  for  the  cylinder  was  not  expressed  in  closed 
form.  Also,  in  both  papers,  the  drag  coefficients  for  the  conical  bodies 
were  calculated  only  for  the  special  case  of  zero  angle  of  attack.  In 
this  paper,  the  drag  coefficients  for  the  sphere  and  cylinder  are  ex¬ 
pressed  in  closed  form,  and  the  calculations  for  the  drag  coefficients 
for  the  cones  are  extended  to  include  angles  of  attack  other  than  zero. 
Two  types  of  molecular  reflection,  diffuse  and  specular,  are  investi¬ 
gated.  Molecular  reflection  which  occurs  in  a  random  direction  without 
relation  to  the  previous  velocity  direction  is  termed  diffuse  reflec¬ 
tion,  Reflection  such  that  the  molecules  leave  the  surface  at  an  angle 
equal  to  the  angle  of  Incidence  is  called  specular  reflection. 

As  previously  stated,  free-molecule  theory  accounts  for  both  the 
mass  motion  of  the  gas  and  the  thermal  motion  of  the  molecules.  Early 
investigations  of  rarefied  gas  flow  were  based  upon  the  assumption  that 
the  velocity  of  the  molecules  was  equal  to  the  translatory  velocity  of 
the  gas,  thus  neglecting  the  thermal  motion  of  the  molecules.  As 
pointed  out  by  Zahm  (reference  6),  expressions  for  the  forces  acting  on 
various  bodies  located  in  a  field  of  such  particles  were  derived  by 
Newton,  who  as  sinned  that  the  molecules  were  either  perfectly  elastic  or 
perfectly  inelastic,  A  perfectly  elastic  molecule  rebounds  after  a 
collision  with  a  surface  in  such  a  manner  that  the  n  rmal  component  of 
velocity  is  reversed.  A  perfectly  inelastic  molecule  loses  the  normal 
component  of  velocity  upon  impact.  In  both  cases,  however,  the  tangen¬ 
tial  component  remains  uncharged.  In  this  paper,  flow  of  elastic  par¬ 
ticles  liaving  only  a  translatory  mass  motion  will  be  termed  "Newtonian 
flow,"  and  flow  of  inelastic  particles  will  be  termed  "inelastic 
Newtonian  flow," 

At  high  flow  speeds  where  molecular  thermal  velocities  become  rela¬ 
tively  unimportant,  the  expressions  derived  from  fre e-molecule— flow 
theory  for  the  forces  on  a  body  for  the  case  of  specular  reflection  re¬ 
duce  to  the  expressions  for  a  body  in  a  Newtonian  flow  field. 

J .".elastic  Newtonian  flow  approximates  the  condition  of  continuum 
hypersonic  flow  in  which  the  shock  wave  lies  very  close  to  the  body. 
Consequently,  impact  forces  of  inelastic  Newtonian  flow  have  been  applied 
to  an  analysis  of  hypersonic  flow  (reference  7).  In  addition  to  the 
impact  forces  of  inelastic  Newtonian  flow,  other  forces  of  a  centripetal 
nature  exist  as  a  result  of  flow  over  curved  surfaces.  An  investigation 
of  the  effects  of  these  centripetal  forces  is  also  presented  in  refer¬ 
ence  7, 
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The  purpose  of  this  paper  is  to  obtain  solutions  for  the  aerody¬ 
namic  coefficients  of  various  bodies  in  a  free-mclecule— flow  field  and 
to  calciilate  the  lift  and  drag  coefficients  for  a  conical  body  at  an 
angle  of  attack,  assuming  two  types  of  molecular  reflection,  specular 
and  diffuse.  A  comparison  is  then  made  of  the  drag  coefficient  for  a 
cone  at  zero  angle  of  attack  calciilated  by  free-molecule  methods,  assum¬ 
ing  a  hypothetical  type  of  molecular  reflection  similar  to  that  of 
Hewton's  inelastic  particles,  with  the  drag  coefficient  tabulated  in 
reference  7  for  inelastic  Newtonian  flow  and  the  drag  coefficient  calcu¬ 
lated  in  reference  8  for  continuum  flow. 


NOTATION 


A  surface  area,  square  feet 


Cy  >  components  of  total  velocity,  feet  per  second 


aerodynamic  coefficient  due  to  impinging  molecules,  dimensionless 

Cj)^  drag  coefficient  due  to  impinging  molecules,  dimensionless 

lift  coefficient  due  to  impinging  molecules,  dimensionless 

Cj.  aerodynamic  coefficient  due  to  diffuse  re-emlsslon  of  the 
molecules,  dimensionless 


drag  coefficient  due  to  diffuse  re-emission  of  the  molecules, 
dimensionless 


Ct  lift  coefficient  due  to  diffuse  re-emission  of  the  molecules, 
^  dimensionless 

Cp  total  drag  coefi’lcient  assuming  diffuse  re-emission 
(Cj)^  +  dimensionless 

total  lift  coefficient  assuming  diffuse  re— emission 
dimensionless 

total  drag  coefficient  assuming  specular  reflection  of  the 
molecules,  dimensionless 


k 
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Gi 

Gi- 

Gs 

lo 

) 

V  ^ 

! 

L 

M 

m 

N 

R 


8r 

U 

Ux  1 
Uy  > 
Uz 


total  lift  coefficient  assimiing  specular  reflection  of  the 
molecules^  dimensionless 

momentim  force  due  to  Impinging  molecule s^  pounds 
normal  diffuse  re-emission  monentum  force ^  pounds 
total  normal  specular  momentum  force,  pounds 
modified  Bessel  function  of  first  kind  and  zero  order 
modified  Bessel  function  of  first  kind  and  first  order 


direction  cosines  in  an  arbitrary  direction,  dimensionless 
(Subscripts  I  and  d  refer  to  lift  and  drag  directions.) 


length,  feet 

Mach  number,  dimensionless 

mass  of  one  molecule,  slugs 

niuaber  of  molecules  per  unit  voIuhb  of  gas 

radius  of  body,  feet 


molecular  speed  ratio  (ratio  of  stream  mass  velocity  to  most 
probable  molecular  speed) 

I  ”) 


dimensionless 


diffuse  re-emission  speed  ratio 
dimensionless 


(4)’ 


masa  velcxjity,  feet  per  second 


components  of  mass  velocity,  feet  per  second 
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Vm  most  protatle  molecular  speed,  feet  per  second 

Yr  most  probable  re-emission  speed,  feet  per  second 

1 

X  I 

i 

;  * 
j  I'  local  Cartesian  coordinates 

z 

a  angle  of  attack  of  body  with  respect  to  free  stream 

3  reciprocal  of  most  probable  molecular  speed 

—  ] ,  seconds  per  foot 
Vm/ 

7  ratio  of  specific  heats,  dimensionless 

6  semivertex  angle  of  cone 

9  angle  of  attack  of  body  element  with  respect  to  free  stream 
p  density  of  gas  stream,  slugs  per  cubic  foot 
^  azimuthal  angle  measured  at  base  of  cone 
X  dimensionless  quantity  defined  by 

X  =  exp  (-s^Ixd^)  +  ^/«'  s  2xd  [1  +  erf  (s  Ixi)] 

X  *  dimensionless  quantity  defined  by 

X*  =  exp  (-s^2^^’^)  —  8  ^xd'  ~  * )] 

dimensionless  quantity  defined  by 

n  =  (2x3Uk  +  V^^y  ^zP^z)  lexpM^U/)  + 

+  “  lx  [1  +  erf  (3Ux)] 
n*  dimensionless  quantity  defined  ty 

fit  =(z^‘pUx*  +  2y'pUy’  +  2z'PUz’)  |exp(-32Ux*2)  _ 

3Ux*  [1  -  erf(BUx*)]|  Ix’  [l  -  erf(3Ux')] 


pUx  [1  +  erf(pUx)] 
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erf (a)  error  function  — 

Superscript 


refers  to  shielded  surface  only 


ANALYSIS 


The  basic  assumptions  involved  in  the  formulation  of  the  expres¬ 
sions  for  the  force  on  a  body  in  a  free-molecule-flow  field  are  the  fol¬ 
lowing:  (l)  the  gas  molecules  have  a  Maxwellian  velocity  distribution 
superimposed  upon  a  uniform  mass  velocity,  and  (2)  collisions  between 
impinging  and  re-emitted  molecules  are  negligible. 

As  a  result  of  the  second  assumption,  the  total  force  can  be  broken 
down  into  two  components ;  one  arising  from  the  bombardment  by  impinging 
molecules  and  the  other  arising  from  the  re-emission  of  the  moleC';les 
from  the  surface. 

The  type  of  re-emlsslon  of  the  moleciaes  from  the  body  surface  can 
be  divided  into  several  categories,  the  more  common  of  these  being  dif¬ 
fuse  and  specular  reflection.  The  phenomenon  of  molecular  reflection  is 
treated  in  greater  detail  in  references  1  and  9,  and  only  a  brief  dis¬ 
cussion  is  given  here.  Diffuse  reflection,  which  is  most  common  physi¬ 
cally,  occurs  in  such  a  manner  that  all  previous  directional  history  is 
erased,  the  actual  direction  of  re-emlsslon  being  controlled  by  the 
Khudsen  cosine  law.  The  nature  of  this  type  of  re-emission  is  such  that 
the  molecules,  upon  leaving  the  surface,  have  a  Maxwellian  distribution 
of  speed  which  depends  upon  the  temperature  of  the  re— emitted  stream. 

For  the  case  of  specular  reflection,  the  molecules  leave  the  surface  in 
a  direction  determined  by  the  angle  of  incidence.  Thus,  the  normal  com¬ 
ponent  of  velocity  is  reversed,  while  the  tangential  component  remains 
unchanged. 

The  method  of  calculation  of  Cj^  and  C^.,  the  incident  and  re- 
emission  aerodynamic  coefficients,  is  outlined  in  reference  4,  and  is 
summarized  in  the  following  paragraphs. 

In  the  Integrals  presented  below,  which  are  taken  from  reference  4, 
the  direction  cosines  between  the  local  coordinate  axes  of  the  body  sur¬ 
face  area  and  the  momentum  force  are  given  bj  2x,  Ij,  and  Iz.  Figure 
1  illustiates  the  coordinate  axes  for  an  element  of  area  used  in  this 
report.  The  positive  x  axis  is  normal  to  the  element  and  is  directed 
into  the  exposed  surface.  The  mass  velocity  vector  makes  an  angle  9 
with  the  positive  direction  of  the  y  axis.  Figures  2  and  3  illustrate 
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the  elements  of  area  chosen  for  each  of  the  bodies*  The  components  of 
mass  velocity  Ux>  Uy,  and  Uz  are  in  the  direction  of  the  local  axes. 

The  first  integral  term  applies  directly  to  the  area  exposed  to  the  mass 
velocity;  the  second  integral  term  applies  to  the  rear  or  shielded  sur¬ 
face  of  the  body.  These  areas  are  taken  into  account  separately  because 
the  front  of  the  body  acts  as  a  shield  preventing  collisions  of  the  rear 
side  with  all  molecules  except  those  having  absolute  velocity  components 
in  the  direction  of  fli^t  with  magnitude  eq.ual  to  or  greater  than  the 
flight  speed  (reference  1).  The  shielded  area  on  the  surface  of  a  cone 
is  shown  in  figure  4*  Other  methods  have  been  devised  to  account  for  the 
shielded  area  and  are  presented  in  references  2  and  3« 


Force  Due  to  Impinging  Molecules 


The  component  of  momentum  in  a  direction  defined  by  Ij,  and 

1 2,  imparted  to  a  differential  plane  area  in  a  fre e-molecule— flow  field 
by  the  impact  of  impinging  molecules  is 


“i  ■;375 


00  nOO  nOO 


Cx(^x^x  +  ^y^y  3 


-oo^O 


{- 


(cx  —  + 


(cj  —  Uy)®  +  (cz  —  Uz)‘ 


dCx  dcy  dCz  - 


Cx(^x®x  ^z^z)  ©xp 


{ 


(cx  —  Ux)^  + 


(Cy  -  Uy)®  + 


(cz 


dCx  dcy  dcz 


dA 


(1) 


or 
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dG.  = 


1 


2  8‘ 


{•/«  [1  +  erf  (pu^)] 


exp  (-P^u/)  [>  +^1£  2^  [1  +  erf  (pu^)] 

I  2 


-/n 


(ix3Ux+W^zP^z){V«  3U3c[1  -  erf  (pu^,)]  -  exp  (-p^U^^)  j 


*/n 


I  Ij  [  1  -  erf  (pUx)] 


dA 


(2) 


vhere  p  =  mN  and  p  =  l/V^  =  s/U. 

The  choice  of  the  direction  cosines  is  related  to  the  component 
of  force  being  computed.  If  the  drag  is  being  calciaatedj  it  is  nec¬ 
essary  to  find  the  direction  angles  between  each  of  the  local  axes 
and  the  direction  of  drag.  Likewise,  direction  angles  are  found  for 
other  force  component  being  investigated. 

In  applying  eq.uation  (2)  to  arbitrary  bodies,  the  momentum  force 
in  a  desired  direction  due  to  the  impinging  molecules  is  first  found 
-for  an  element  of .area  of  the  body  under  Investigation.  When  this 
relationship  is  used,  an  Integration  over  the  surface  area  of  the  body 
gives  the  desired  momentum  force  on  the  entire  body  due  to  Impinging 
molecules.  The  incident  aerodynamic  coefficient  referred  to  a  char¬ 
acteristic  body  area  is  then  found. 


Force  Due  to  Diffuse  Molecular  Deflection 


For  the  case  of  diffuse  molecular  reflection,  the  normal 
re-emission  force  on  the  front  side  of  a  differential  area  inclined 
to  the  direction  of  flow  is  given  by 


dG^  =  ^  _2L_  clA 


2  2ss., 


(3) 


and  the  normal  re— emission  force  on  the  rear  side  is  given  by 


dG,. 

^  2  2ssj. 


(h) 
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The  difference  obtained  by  subtraction  of  equation  (4)  from  equation  (3) 
is  the  total  re-emission  force .  (See  appendix  A,  reference  4. ) 

The  value  of  Sj.  is  a  complicated  function  of  the  surface  condi¬ 
tion  of  the  body;  however,  in  this  paper,  Sj*  is  assumed  equal  to  s 
for  all  calculations.  At  low  flow  speeds,  this  assumption  is  very 
nearly  correct.  At  hi^  flow  speeds,  however,  the  value  of  s^  is  less 
than  that  of  s. 

The  resolution  of  the  normal  force  of  equations  (3)  and  (4)  into 
the  lift  and  drag  components  followed  by  integration  over  the  body  sur¬ 
face  gives  the  re-emission  lift  and  drag  acting  on  the  body.  After 
calculation  of  the  lift  and  drag,  the  aerodynamic  coefficients  for  dif¬ 
fuse  molecular  reflection  can  be  found. 


Force  Due  to  Specular  Moleculsir  Reflection 


For  the  case  of  specular  reflection,  the  molecules  leave  the  sur¬ 
face  in  a  direction  determined  by  the  angle  of  Incidence.  Since  the 
reflection  angle  is  equal  to  the  angle  of  incidence,  the  normal  component 
of  velocity  is  reversed,  while  the  tangential  component  remains  un¬ 
changed.  Therefore,  the  momentum  change  is  in  a  direction  normal  to  the 
bombarded  surface,  and  the  total  momentum  force  consists  of  two  equal 
components:  one  due  to  impingement  of  the  molecules,  the  other  due  to 
the  reflection  of  the  molecules  from  the  surface. 


The  incident  momentum  force  on  a  differential  area  in  a  direction 
normal  to  the  surface  is  found  from  equation  (2)  by  the  substitution  of 
the  following  direction  cosines: 

=  cos  0=1 

2y  =  cos  2=0  (5) 

=  cos  —  =  0 
z  2 

The  total  normal  momentum  force  for  the  case  of  specular  reflection  is 
twice  as  great  as  the  incident  normal  force,  and  it  is  given  by 


—  pU  2 


^/^^[l  +  erf(pUj)]  2}  pUxexp(H3®Ux®)|  - 

I  [1  -  erfOUx)]  “  PUxexp (-3^Ux2 )|  dA  (6) 
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Integration  of  equation  (6)  over  the  surface  of  the  body  after  the  nor— 
mal  force  has  "been  resolved  into  the  lift  and  drag  components  yields  the 
total  lift  and  drag  on  the  body  due  to  specular  reflection. 


Force  Due  to  Molecular  Reflection  Analogous  to  Hypersonic  Flow 


A  third  type  of  molecular  reflection  is  postulated  here  in  order  to 
approximate  the  characteristics  of  hypersonic  flow  over  a  body.  This 
type  of  reflection  is  much  like  that  of  Newton*s  inelastic  particles^ 
the  only  difference  being  that  thermal  velocities  are  considered.  It  is 
assumed  that  the  molecules  flow  along  the  body  after  collision  with  the 
surface;  thus^  the  normal  component  of  velocity  is  destroyed  while  the 
tangential  component  remains  unchanged.  The  values  of  the  aerodynamic 
coefficients  thus  obtained  are  one-half  the  values  calculated  for  specu¬ 
lar  reflection.  This  type  of  molecular  reflection  is  extremely  unlikely 
and  probably  physically  impossible;  however^  the  Justification  for  this 
assumption  lies  in  the  fact  that  the  conditions  of  hypersonic  continuum 
flow  over  a  body  are  approximated.  In  continuum  flow^  the  stream  fol¬ 
lows  the  surface  of  the  body.  As  the  flow  speed  increases,  the  zone  of 
influence  of  the  body  on  the  gas  stream  approaches  the  surface  of  the 
body  until  the  shock  wave  lies  nearly  on  the  body  surface#  On  the  basis 
of  the  above  assumption,  free-molecule  theory  is  compared  with  continuum 
theory. 

The  expressions  for  the  aerodynamic  coefficients  of  a  flat  plate, 
cylinder,  sphere,  and  cone  together  with  the  details  of  the  development 
are  presented  in  the  appendix. 


DISCUSSION  AND  RESULTS 
Flat  Plate 


The  aerodynamic  coefficients  for  a  flat  plate  Inclined  at  various 
angles  of  attack  to  the  direction  of  mass  flow  were  calculated  for 
values  of  molecular  speed  ratio  varying  from  0  to  20  and  are  presented 
in  figures  5,  6,  and  7. 

The  limiting  values,  as  the  molecular  speed  ratio  approaches  infin¬ 
ity,  of  the  aerodynamic  coefficients  for  a  flat  plate  for  the  case  of 
specular  reflection  reduce  to  the  classical  expressions  for  Newtonian 
flow,  mentioned  by  Zahm  (reference  6) 

Cd^  =  h  sin®  a 

=  h  sln^  a  cos  a 
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A  ccmpaxison  of  the  drag  coefficlentg  for  diffuse  spec^llar 
reflection  indicates  that  the  value  for  diffuse  reflection  has  a  larger 
value  at  low  angles  of  attack;  at  the  hi^  angles  of  attack,  the  spec¬ 
ular  drag  coefficient  is  greater.  This  is  explained  hy  the  fact  that, 
at  small  angles  of  attack,  the  drag  for  the  case  of  specular  reflection 
is  due  to  the  small  normal  component  of  velocity;  whereas  the  drag  for 
the  case  of  diffuse  reflection  has  a  value  which  results  from  the  loss 
of  both  a  small  normal  velocity  component  and  a  large  tangential  velo¬ 
city.  At  large  angles  of  attack,  the  re— emission  velocity  arising  from 
diffuse  reflection  is  less  than  the  re-emission  velocity  due  to  speciilar 
reflection,  while  the  momentum  effects  of  the  incident  molecules  are 
nearly  the  same  for  both  cases.  Thus,  the  total  drag  coefficient  due  to 
specular  reflection  is  greater.  This  last  point  is  made  clear  by  consid¬ 
ering  a  flat  plate  at  a  90°  angle  of  attack. 

The  lift-drag  ratios  for  the  plate  are  shown  in  figure  6.  The  ex¬ 
pression  for  the  lift— drag  ratio  for  specular  reflection  is  cot  a. 
Obviously,  these  values  have  a  great  range.  The  aerodynamic  coefficients 
for  double-^edge  wings  can  easily  be  found  by  modifying  the  coefficients 
for  the  flat  plate.  The  low  values  for  diffuse  reflection  are  notable. 


Cylinder 


Values  of  the  total  drag  coefficient  for  the  case  of  diffuse  re¬ 
flection  were  obtained  from  reference  4  and  are  presented  in  figure  8. 
An  analytical  expression  for  Sj.  employed  in  the  calculations  of  ref¬ 
erence  4  was  found  by  means  of  an  energy  and  mass  balance  over  the  sur¬ 
face  of  the  cylinder  taking  into  account  conditions  which  exist  at  the 
surface  of  the  cylinder.  Details  of  the  calculations  can  be  found  in 
appendix  A,  reference  4. 

A  recent  experimental  investigation  of  drag  (reference  4)  has  sup¬ 
plied  the  drag-coefficient  data  plotted  in  figure  8.  Theoretical 
curves  of  the  drag  coefficient  for  both  diffuse  and  specular  reflection 
are  plotted  for  comparison  with  the  experimental  data.  A  con^jlete  dis¬ 
cussion  of  the  comparison  of  the  experimental  data  and  the  curve  repre¬ 
senting  diffuse  reflection  can  be  found  in  reference  4.  In  addition,  a 
curve  of  the  drag  coefficient  for  the  case  of  diffuse  reflection,  which 
includes  the  assumption  that  the  molecular  speed  after  impact  is  the 
same  as  that  before  impact,  is  shown  for  comparison  with  the  drag  co¬ 
efficient  obtained  from  reference  4,  The  curves  agree  well  at  the  low 
speeds,  but  diverge  as  the  molecular  speed  ratio  increases. 
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Sphere 


The  results  of  the  calculations  for  the  drag  coefficient  for  a 
sphere  are  presented  in  figure  9*  The  drag  coefficient  for  the  case  of 
diffuse  reflection  is  coinpared  with  the  drag  coefficient  for  the  case  of 
specular  reflection.  The  value  of  the  drag  coefficient  for  diffuse 
reflection  lies  above  the  value  for  specular  reflection  but  both  curves 
asymptotically  approach  the  limiting  value  of  2  at  large  values  of  s. 

The  investigation  into  the  case  of  specular  reflection  for  a 
sphere  yields  an  Interesting  result.  Equations  (AI5)  and  (AI9)  indi¬ 
cate  that  the  total  drag  due  to  the  specular  type  of  reflection  has  the 
same  value  as  the  drag  caused  by  molecules  striking  the  surface  and 
coming  to  rest;  that  is,  the  integrated  effect  of  the  re-emlsslon  drag 
arising  from  specular  reflection  of  the  molecules  from  the  surface  is 
zero.  It  should  be  noted  that  the  sphere  la  unique  in  that  the  re- 
emission  momentum  force  of  the  molecules  reflected  in  a  forvard  direc¬ 
tion  is  equal  and  opposite  to  the  momentum  force  of  the  molecules  re¬ 
flected  in  the  direction  of  the  mass  flov. 


Cones 


Calculations  for  the  drag  and  lift  coefficients  of  the  various 
cones,  assuming  three  types  of  molecular  reflection,  are  shown  in 
figures  10,  11,  12,  I3,  and  l4. 

Results  considering  diffuse  reflection.—  Figure  10  gives  the  drag 
coefficient,  Cp  =  Cp^  +  Cp  ,  for  the  case  of  diffuse  molecular  reflec¬ 
tion.  It  is  interesting  to  note  that  the  drag  coefficients  for  the 
sharp  cones  are  higher  than  those  for  the  blunt  cones.  The  impinging 
molecules  lose  both  the  tangential  and  normal  components  of  velocity. 

The  tangential  or  friction  force  on  the  greater  exposed  surface  area  of 
the  sharp  cones  is  greater  than  the  tangential  force  on  the  blunt  cones. 
Hence,  the  drag  coefficients  are  large  for  the  sharp  cones  since  the 
reference  area  is  the  same  base  area  for  all  cones. 

The  lift  coefficient,  Cp  =  Cp^  +  Cp^,  is  plotted  in  figure  11. 

The  lift  coefficient  for  blunt  cones  at  large  angles  of  attack  de¬ 
creases  as  molecular  speed  ratio  becomes  small.  This  is  due  to  the  col¬ 
lision  of  the  top  or  shielded  surface  of  the  cone  vlth  molecules  posses¬ 
sing  an  absolute  velocity  in  the  direction  of  flight  with  magnitude 
equal  to  or  greater  than  the  flight  speed.  The  momentum  force  due  to 
these  molecules  is  then  in  a  direction  opposite  to  the  direction  of 
lift. 
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At  low  values  of  molecular  speed  ratio,  in  the  range  where  the 
thermal  velocity  is  equal  to  or  greater  than  the  mass  velocity,  the 
contribution  of  the  thermal  velocity  to  the  absolute  velocity  becomes 
a  determining  factor  in  the  calculation  of  the  aerodynamic  coeffi¬ 
cients;  whereas,  at  high  values  of  molecular  speed  ratio,  the  contri¬ 
bution  of  the  thermal  velocity  is  negligible. 

Ee suits  considering  specular  reflection.-  Figures  12  and  I3  give 
the  drag  and  lift  coefficients  for  cones  for  the  case  of  specular  re¬ 
flection.  At  small  angles  of  attack,  the  drag  coefficients  of  the  blunt 
cones  have  higher  values  than  those  of  the  sharp  cones.  The  reverse  is 
true  at  large  angles  of  attack.  A  very  sharp  cone  at  zero  angle  of 
attack  has  very  little  drag  because  the  change  in  momentum  of  the  mole¬ 
cules  has  only  a  very  small  component  in  the  direction  of  drag.  This 
component,  acting  on  a  blunt  cone,  is  much  greater.  At  a  high  angle  of 
attack,  the  surface  area  exposed  to  the  flow  by  the  sharp  cones  exceeds 
■the  exposed  surface  area  of  the  blunt  cones,  for  equal  base  areas,  and 
hence  the  drag  force  is  greater  for  the  sharp  cones.  Since  the  drag 
coefficients  are  referred  to  the  sane  base  area  rather  than  a  projection 
of  the  exposed  area,  they  must  vary  in  direct  proportion  to  the  drag 
force . 


Eesults  considering  an  approximation  of  continuum  theory.—  The 
values  of  the  lift  and  drag  coefficients  obtained  by  means  of  the  hypo¬ 
thetical  type  of  re— emission  postulated  in  the  analysis  are  half  the 
values  obtained  for  the  case  of  purely  specular  reflection  in  which  the 
nonnal  component  of  velocity  is  reversed  as  a  molecule  leaves  the  sur¬ 
face  . 


Values  of  drag  coefficients  for  cones  in  continuum  flow  have  been 
calculated  by  Kopal  (reference  8).  Likewise,  reference  7  describes  an 
investigation  of  continuum  flow  at  hypersonic  speed.  In  figure  l4,  the 
drag  coefficient  of  a  30°  semivertex  angle  cone,  at  zero  angle  of 
attack,  calculated  for  the  above  hjrpothetical  case  of  molecular  re¬ 
flection,  is  compared  with  the  drag  coefficient  calculated  in  references 
7  and  8. 

Experimental  values  for  the  drag  coefficient  of  a  30°  semivertex 
cone-cylinder  fired  in  the  Ames  supersonic  free— flight  tunnel  (reference 
10)  are  also  shown  in  figure  l4  for  comparison  with  the  analytic 
values.  The  experimental  values  include  head  drag  on  the  cone,  friction 
drag  on  the  cylindrical  afterbody,  and  base  drag.  The  free-molecule 
values  take  into  account  head  drag  on  a  30°  semivertex  angle  cone  only. 
At  high  values  of  flow  speed,  the  drag  coefficient  values  obtained  from 
free-molecule  theory  approach  the  values  for  inelas"'  ic  Newt.onlan  flow. 
Although  the  analysis  of  reference  7  foi"  flow  past  a  cone  takes  into 
account  centrifugal  forces  arising  from  flow  over  curved  surfaces,  the 
expression  for  the  drag  coefficient  reduces  to  that  for  inelastic  New¬ 
tonian  flow  at  zero  angle  of  attack  where  the  centrifugal  forces  are 
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zero.  The  drag  coefficient  due  to  inelastic  Newtonian  flow  is  repre¬ 
sented  in  figure  l4  by  a  straight  line,  since  it  Is  independent  of  flow 
speed. 

In  figure  l4,  the  agreement  between  free-molecule— flow  theory  anri 
the  experimental  data  for  a  continuum  flow  at  low  values  of  molec^llar 
speed  ratio  is  entirely  fortuitous.  The  increase  in  the  value  of  the 
drag  coefficient  calculated  from  free— molecule— flow  theory  is  due  to  the 
fact  that  the  drag  force  is  proportional  to  the  first  power  of  the  velo¬ 
city  at  low  speeds;  whereas  the  increase  in  the  value  of  the  experimental 
drag  coefficient  is  the  characteristic  rise  which  occurs  in  the  tran¬ 
sonic  range  of  stream  speeds  in  continuum  flow. 


CONCLUDING  EEMAEKS 


The  lift-drag  ratios  of  a  flat  plate  in  free-molecule  flow  are 
very  small  for  the  case  of  diffuse  molecular  reflection.  The  lift- 
drag  ratios  decrease  as  the  molecular  speed  ratio  is  increased.  At 
high  speeds,  the  lift-drag  ratios  for  the  flat  plate  are  more  favorable 
at  snail  angles  of  attack  than  at  large  angles  of  attack. 

For  the  case  of  diffuse  molecular  reflection,  the  lift  and  drag 
coefficients  of  all  the  bodies  investigated  approach  constant  limiting 
values  at  hi^  speeds.  These  limiting  values  depend  upon  the  body 
geometry  and  the  angle  of  attack.  At  low  speeds,  the  drag  coefficients 
for  all  the  bodies  approach  an  infinite  value;  this  is  due  to  the  fact 
that  the  drag  force  is  proportional  to  the  first  power  of  the  velocity 
at  low  speeds. 

For  the  case  of  specular  reflection,  the  aerodynamic  coefficients 
of  all  the  bodies  investigated  likewise  approach  constant  limiting 
values  at  hi^  speeds.  These  limiting  values  correspond  to  those  cal¬ 
culated  by  means  of  inelastic  Newtonian  flow  theory. 


Ames  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Moffett  Field,  Calif.,  May  9,  1951. 
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APPENDIX 

DEEIYATION  OF  THE  AERODYNAMIC  COEFFICIENTS  FOR 
BODIES  IN  A  FREE-MOLECULE-^OW  FIELD 

The  details  of  the  development  for  the  lift  and  drag  coefficients 
for  the  various  bodies  are  given  below; 

Flat  Plate 


The  combined  form  of  equation  (2)  ia 


^  Z  Zy  pTJy+Z  z  pTJz 


) 


pUx  erf(pUx)+exp(-P^Ux^) 


+ 


(Al) 


Since  the  faces  of  a  flat  plate  are  parallel^  the  same  direction 
cosines  apply  to  hoth  the  front  and  rear  surfaces  and  the  momentum  of 
the  molecules  striking  hoth  faces  can  he  collectively  taken  into  account, 
and  equation  2  written  in  the  comhined  form  is  appllcahle-  This  is  not 
true  for  an  element  of  area  chosen  for  the  cone,  however,  and  the  com¬ 
bined  form  of  equation  (2)  cannot  he  used  for  that  case.  For  the  case 
of  the  flat  plate,  equation  (Al)  gives  the  incident  momentum  force  on 
hoth  sides  of  the  plate,  and  the  aerodynamic  coefficients  can  he  found 
without  further  Integration. 


The  direction  cosines  and  velocity  components  are  the  following: 


Zxd  =  sin  a 

=  —COS  Ot/ 

Z  zd  ~  ^ 


ZxZ  =  cos  a 
Ijl  =  sin  a 
^zZ  =  0 


Ux  =  U  sin  a 
XJy  =  -U  COS  a 
Uz  =  0 


(A2) 


3U  =  3  J 


Substitution  of  the  relations  of  equation  (A2)  into  equation  (Al) 
gives  the  following  coefficients,  referred  to  the  area  of  one  side  of 
the  plate. 
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2  2 

exp  (-S  sin  a)+2  sin  a 


erf  (s  sin  a)  (A3) 


cos  a 


erf  (s  sin  a) 


(A4) 


The  normal  force  due  to  diffuse  molecular  re— emission  from  both 

front  and  rear  surfaces  is  given  by  the  difference  between  equations  (3) 
and  (4).  1 


(iG-27  ““  d.CT];7  * 


pu^  1 

2  2ssr 


(X  -X')  dA 


(A5) 


Substitution  of  U  sin  a 
reduces  eq[uation  (A5)  to 


dGj^  —  dGp  * 


for  Ux  and  s  for  PU  in 


PU2  sin  a 


dA 


X 


and  X« 


(A6) 


The  re—emission  drag  force  on  an  element  of  area  is 

pTj^  sin^a 
_  _ 


dA 


and 


'Dr 


sin  a 


Sr 


(A?) 


The  total  drag  coefficient  is  then 
=  — ^exp  (-s^  Sin^a)  +  2  sin  a  (1  +  erf  (s  sin  a)+ 


'D 


\/k  b 


The  lift  due  to  molecular  re-^mission  on  an  element  of  area  is 

sin  g  cos  a 

2  St* 


dA 


and 


% 


y^  sin  a  cos  a 


(A8) 


The  total  lift  coefficient  is  then 


cos  a _ ^  ^  a/jT  sin  a  cos  a 

- -  erf  (s  sin  a)  +  -  ■  - .  . 
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For  the  case  of  specular  reflection,  the  total  momentum  force  in  a 
direction  normal  to  the  surface  resulting  from  incident  and  reflected 
molecules  on  both  sides  of  a  flat  p)late  is 


dGa  — 


plje  2 


2  JlTs  L 

from  eq,uation  (6). 


2PUxexp(  ^  erfOUx)  dA 


The  drag  and  lift  coefficients  are  then  found  after  resolving  the 
normal  force  as  given  in  equation  (A9)  into  drag  and  lift  components. 

h 

CDg  =  sin  a  - — •  sin  a  exp(-s2sin2a)  +  1 


^sin^a  +  ^  erf(s  ^in  a) 


(AlO) 


=  cos  a 


sin  a  exp(— s^sin^a)  + 


s 

4  ^sin^a  +  erf(s  sin  a) 


(All) 


Cylinder 


The  case  of  diffuse  molecular  reflection  on  a  transverse  cylinder 
has  been  investigated  in  reference  4,  and  the  total  drag  coefficient 
based  on  the  projected  area  of  the  cylinder  is 


(l+2s2) 


.(?) 


The  specular  case  only  remains  to  be  considered  here, 
of  area  is  chosen  as  shown  in  figure  2,  thus 

dA  =  L  R  d  0 


The  element 


^xd  =  ® 

lyd  =  -COS  0 


Uj.  =  U  sin  0 
Uy  =  -U  COS  0 


l.A  =  0 


U,  =  0 
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Each  element  of  area  Is  considered  to  have  a  rear  or  shielded  area 
and  the  normal  force  on  the  cylinder  is  found  from  equation  (A9). 


dGs 


pu^  4lr  r  /  ^  .  X 

=  - - — -  I  2a  sin  6  exp(-s2sin2e)  + 


2  x/rt  ^s^sin^e  +  erf(8  sin  0)J  d© 


The  drag  force  is 

pTJ^  4tR 


2  V 


sin  0 


4lr  r 

2  ^s^sin^©  +  erf(8  sin  6)jd0 


2s  sin  0  exp(-e^sin^0)  + 


(A12) 


The  integrals  are 


8in20  exp(-B^sin20)  40  =  ^  exp  ^ 


«/ 


J'  sin  0  erf  (s  sin  0)  d0  =  exp  (^1”^  ■*■  (1”^ 


*/2 


sin®0  erf  (s  sin  0)  d0  =  — ^ -  exp  ^^o 

(f)] 

and  the  drag  coefficient  due  to  specular  reflection^  referred  to  the 
frontal  projected  area^  is 


ft  _  2  Vjt  / 

%  "  3b  ^  2  ^ 


(3+232)  +  (i+2s2)  (_) 


(A13) 


Sphere 


As  shown  in  figure  2,  the  element  of  area  for  the  sphere  is  given 


t>7 


dA  =  2rtR2  cos  0  d0 


WACA  TN  2k23 


19 


Equation  (Al)  for  tlie  incident  momentum  drag  force  is 

P'0®  /  2«  2  cos  e  ,  2  P  . 

dGi  =  (itR  )  -  23  eip(-e  sin  0)  + 

-/IT  B®  ^ 


yit  sin  0  (2b2  +  l)  erf(s  sin  0)  d0 


The  integrals  are 


(Al4) 


B  cos  0  exp(— s®sin®0)  d0  =  — ^  erf(s) 


cos  0  sin  0  erffs  sin  0)  d0  =  ■  ■-  ‘ 

'  ^  2^s  2 


erf (s ) 


and  the  incident  drag  .coefficient,  referred  to  the  frontal  projected 
area,  is 


(A15) 


The  drag  force  due  to  molecular  re-emission,  from  equation  (A6),  is 


po®  r  2^ 


Bin^O  cos  0  d0 


O  ^3? 


(A16) 


Integration  of  equation  (Al6)  yields  the  re— emission  drag  coefficient 

V  /  A  -1 


'Or  3sr 


(A17) 


The  total  drag  Is  given,  for  diffuse  reflection,  hy 

^  2  exp(-e^)  ,  1  \  1  1  \  \  2^/l^ 

Cj) - -  (1  +  -— )  +  2(1  +  —  — j-)erf(s)  +  — — 

,/it  S  28^  3^  4s^  3sr 


The  drag  force  due  to  specular  reflection  on  the  sphere  Is  given  hy 


4  008  0  Sin  0  23  3^^  Q  exp(-«23in20)  + 

3®  L 


2  ^^s2sin®0  +  ^  ^  erf(s  sin  0)  de 


(Al8) 
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The  integrals  are 


6  3in^9  exp(— s^ain^O)  d0  =  erf  (a) - ^ 

43^  23^ 

0  sin  0  erf  (3  sin  0)  d0  =  erfCs) 

2  Vn  s  \  ka^  / 


n 

I 


1t/2 


COS  0  3in®0  erf(8  sin  0)  (10  =  '—■  eip(-6^)  +  erf(s) 

S^/lFs®  163* 


and  the  drag  coefficient  due  to  specular  reflection,  referred  to  the 
frontal  projected  area,  is 


23^+1 
J~rc  3^ 


exp (-6^) 


Cone 


(A19) 


Special  attention  must.-be  given  to  the  conical  body.  Two  separate 
cases  must  be  investigated,  depending  on  the  relationship  of  the  angle 
of  attack  of  the  cone  and  the  semi  vert ex  angle  of  the  cone .  At  an  angle 
of  attack  less  than  or  equal  to  the  semi  vertex  angle,  the  entire  surface 
area  is  exposed  to  the  mass  velocity,  and  the  total  momentum  force  is 
evaluated' by  a  single  integration.  For  those  instances  where  the  angle 
of  attack  exceeds  the  semi  vertex  angle,  a  portion  of  the  surface  area  is 
shielded  from  the  mass  velocity,  and  two  separate  integrations,  one  for 
the  frontal  exposed  area  and  one  for  the  shielded  area,  are  necessary  to 
compute  the  momentum  force,  i Aerodynamic  coefficients  are  then  found 
referred  to  the  base  area  of; the  cone.  The  element  of  area  chosen  for 
constant  angle  of  attack  as  shown  in  figure  3  is  given  by 


dA  = 


^  dcp 
C03  6  2 


Case  I  —  Angle  of  attack  less  than  or  equal  to  the  aeml vertex 
angle  The  aerodynamic  coefficient  due  to  impinging  molecules  ref^erred 
to  The  base  area  of  the  cone  la  given  by  (from  the  first  term  of 
equation  (2)) 
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(A20) 


where 


=  (^lxPUx+lyPUy+22pU2^|exp(-p^^2)  +  erf(3U3;.)]  |  + 


yit  lx 


[1  +  erf(pUy)] 


(A21) 


Suhstitution  of  the  appropriate  direction  cosines  into  equation 
(A20)  yields  the  corresponding  lift  or  drag  coefficient. 


For  drag  the  direction  cosines  are 


^xd  ”  ^  ^  ““  ^  ^  ^ 

=  “  cos  a  sin  5  —  sin  a  sin  &  cos  cp 

^zd  ^  ^ 


(A22) 


and  for  lift  the  direction  cosines  are 


^xZ  ”  ”  ^  ^  ^  ^  ^  ^ 

Zy2  =  sin  a  cos  b  —  cos  a  sin  6  cos  cp 


(A23) 


The  expressions  for  the  several  velocity  components  are  the 
following: 


(  (A24) 


=  U(cos  a  sin  6  —  sin  a  cos  6  cos  cp) 
Uy  =  — U(cos  a  cos  D  +  sin  a  sin  6  cos  cp) 
Ug  =  U  sin  a  sin  cp 


At  zero  angle  of  attack,  equation  (A20)  reduces  to  a  readily 
integrahle  expression  and  the  aerodynamic  coefficients  for  this  case 
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and 


3  exp(-e2sln^6)  +  sin  6 


(3^  +  i  ) 


1  +  erf(s  sin  6) 


s^sln  6 


(A25) 


Clj  =  0 

The  lift  and  drag  coefficients  due  to  diffuse  molecular  reflection 
are  found,  from  equation  (3),  to  he 


2it8Sj,sln  5  Zxd  ^<3cp 


(A26) 


and 


2jtss 


1  f 

rsln  6Jq 


(A27) 


At  zero  angle  of  attack. 


Cd  = 


r  2ss- 


(A28) 


and 


Cl^  =  0 


The  following  relations  hold  for  the  total  lift  and  drag  coeffi¬ 
cients  for  the  case  of  specular  reflection: 

It 

*^1*3  ~  “1/2  2  ^  f  ^xd  exp(-p^x^)  + 

^  rtS/^s^sin  6j„  I 

«/ir 


1  +  erf(pUx)  I  (  +  ^J  JdT 


(A29) 


and 
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For  zero  angle  of  attack 

Cjj  =  ^  ^  exp(— s^sln^b)  +  2  faln^  d  +  [l  +  erf  (a  sin  6)] 

2  a/«  a  V  2s^y 

(A31) 

and 

Case  II  —  Angle  of  attack  greater  than  the  aemlTertea:  angle,—  In 
thia  caae,  part  of  the  aurface  area  ia  ahielded  from  the  main  flow. 

Thia  ahielded  area  ia  homded  hy  atralght— line  elementa  of  the  cone 
and,  at  theae  hovmdariea,  the  flow  ia  tangent  to  the  aiur^ace  of  the 
cone;  that  la,  the  component  of  velocity  in  a  direction  normal  to  the 
aurface  ia  zero. 

From  equation  (A2i|-) 

Uj.  =  U(co3  a  ain  &  —  aln  a  coa  5  coa  <pi)  =  0 

at  cp  =  cpi  (A32) 

and 


tan  6 

cos  - 

tan  a 


where  cpi  denotes  the  hoimdary  line  between  the  exposed  and  shielded 
surface  areas  and  is  measured  around  the  cone  from  the  vertical  as 
shown  in  figure  4. 


The  total  Incident  momentum  for  the  entire  surface  area  is 


2k 
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Gi  = 


pU^  tan  & 


,cpi 


2  3“C03  6 


(Z3;.‘PUx‘+ly‘PUy«+l2*PU2')  |  6  X?  (  .p  ^  ^  )  - 


pUx*  [1  -  erf(pUx')]l'  “ ^x’  [l  -  erf  (pUx' )  ] 


dcp  + 


'cpi 


(lxPUx+iyPUy+l2pU2  )  |exp(-p^^^)  +^/rt  pu^;.  [1  +  erf(pUx)]  |  + 


X  [1  +  erf(pUx)] 


dcp 


(A33) 


vhere 


7.  *  -  — ? 

= 

-Ux 

V  "  V 

Uy*  = 

% 

► 

(A34) 

=  -^z 

Uz’  = 

-Uz 

4 

A  change  In  the  coordinate  axis  of  the  element  of  area  on  the 
shielded  surface  of  the  cone  is  the  reason  for  the  use  of  the  primed 
quantities.  The  positive  direction  of  the  x  axis  for  the  exposed  area 
is  into  the  surface;  whereas  the  positive  direction  of  the  x  axis  for 
the  shielded  area  is  out  of  the  surface.  Since  the  direction  cosines 
have  "been  derived  for  the  exposed  area,  reversal  of  the  x  axis  over  the 
shielded  area  results  in  the  negative  quantities.  For  the  previous 
bodies  considered,  this  procedure  was  unnecessary  since  both  front  and 
rear  areas  were  jointly  taken  into  accoimt. 


The  solution  of  equation  (A33)j  for  the  incident  aerodynamic 
coefficient  referred  to  the  base  area,  yields 


C 


i 


_ 1 _ 

^3/2g2gin  g 


(A35) 


where  H  and  are  the  integrands  of  equation  (A33)-  Either  the 

lift  or  drag  coefficients  can  be  calculated  by  substituting  the  appro¬ 
priate  direction  cosines. 


The  re-emission  drag  and  lift  coefficients  for  diffuse  reflection 

are 


1 

2TCssj.sin  6 


Zxd’  x»dcp  + 


(A36) 


'o 
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X’dCp  +r  ZxlXd<P 


2jtssj.sin  & 


The  drag  coefficient  for  specular  reflection  is  expressed  as 


(A37) 


-iW'  exp(-p%'^)  - 


v5r  [1  -  erfOU^’ )]  +  I  dtp  + 

ixd|pUx  exp(-32Ux2)  [1  +  erfOUx)]  (p^x^  +  |^}dcp  (A38) 


The  total  lift  coefficient,  ^^3^  found  by  substitution  of  2^2 
^*037  2xd*  Th©  values  of  the  above  integrals  were  found  by  means  of 
numerical  methods. 
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Figure  5.-  Flat -plate  drag  coefficient,  diffuse  reflection. 


Ratio  of  lift  and  drag  coefficients^ 
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Figure  6.- Fiat -plate  lift- drag  ratio,  diffuse  reflection. 
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Figure  7-  Fiat- plate  drag  coefficient,  specular  reflection. 
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Figures.-  Sphere  drag  coefficient^  diffuse  and 

specular  reflection. 
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Molecular  speed  ratio,  s 


Figure  14.-  Cone  drag  coefficient  comparison  of  free 
molecule  and  continuum -flow  theory;  a*  01 
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